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Abstract 

Background: Somatic mosaicism of copy number variants (CNVs) in human body organs and de novo CNV event 
in monozygotic twins suggest that de novo CNVs can occur during mitotic recombination. These de novo CNV 
events are important for understanding genetic background of evolution and diverse phenotypes. In this study, we 
explored de novo CNV event in cloned dogs with identical genetic background. 

Results: We analyzed CNVs in seven cloned dogs using the nuclear donor genome as reference by array-CGH, and 
identified five de novo CNVs in two of the seven clones. Genomic qPCR, dye-swap array-CGH analysis and B-allele 
profile analysis were used for their validation. Two larger de novo CNVs (5.2 Mb and 338 Kb) on chromosomes X 
and 19 in clone-3 were consistently validated by all three experiments. The other three smaller CNVs (sized from 
36.1 to76.4 Kb) on chromosomes 2, 15 and 32 in clone-3 and clone-6 were verified by at least one of the three 
validations. In addition to the de novo CNVs, we identified a 37 Mb-sized copy neutral de novo loss of heterozygosity 
event on chromosome 2 in clone-6. 

Conclusions: To our knowledge, this is the first report of de novo CNVs in the cloned dogs which were generated by 
somatic cell nuclear transfer technology. To study de novo genetic events in cloned animals can help understand 
formation mechanisms of genetic variants and their biological implications. 

Keywords: Cloned dog, De novo copy number variation, Mosaicism, Loss of heterozygosity, Somatic cell nuclear transfer 



Background 

Together with single nucleotide polymorphisms (SNPs), 
DNA structural variations generally termed copy num- 
ber variations (CNVs) are the major part of genomic var- 
iations in diverse animals including dogs [1,2]. CNVs are 
thought to contribute to phenotypic diversities including 
disease susceptibility and provide the important sub- 
strate for evolution [3,4]. Although it has been widely 
believed that all the cells in one individual are genetically 
identical, this dogma is changing due to our understand- 
ing of genomic variants. Indeed, recent studies have 
demonstrated the somatic mosaicism of CNVs in body 
organs [5,6], which suggests that de novo CNV events 
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may commonly occur during mitotic recombination [7]. 
This possibility was clearly demonstrated by Bruder et 
als report identifying de novo somatic CNV events in 
monozygotic (MZ) twins [8]. They showed that the de 
novo post-twinning CNV frequency could be as high as 
5% on a per-individual basis or 10% per twinning event. 
Also, Breckpot et al. reported CNVs in 1 out of 6 pheno- 
typically discordant MZ twins [9]. These results sug- 
gest that CNV analysis in phenotypically discordant MZ 
twins can be useful to identify genetic loci associated 
with various traits. However, the de novo rate of somatic 
CNV mosaicism events is still largely unknown. For the 
case of cloned dogs, there has been no report about the 
de novo somatic CNV mosaicism event itself. 

Dogs (Canis familiaris) have been used as working an- 
imals due to their superior sniffing ability to any other 
species or machine. For example, a cancer-sniffing dog 
can detect the existence of colorectal cancer from stool 
samples with 97% sensitivity and 99% specificity [10]. 
However, up to 70% of dogs that were bred from 
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working dogs are not suitable for practical use because 
they lack trainability [11]. To reduce the tremendous 
cost, effort, and time to train one working dog, dog- 
cloning technology has been applied for the propagation 
of elite working dogs [12]. 

In this study, we aimed to identify de novo somatic 
CNV events which may exist in cloned dogs with identi- 
cal genetic backgrounds like in human MZ twins using 
microarray-based comparative genomic hybridization 
(array-CGH). Array-CGH is one of the most popular 
and useful tools to analyze genome-wide CNV profiles 
and a number of studies have identified canine CNVs 
using this technique [13-17]. We validated the CNVs by 
diverse analyses such as genomic qPCR, dye-swap array- 
CGH, and B-allele profile analysis. We examined the 
genome-wide de novo somatic CNV events for the seven 
cloned dogs generated in our previous study [12] and 
identified five de novo CNV events. 

Results 

Validation of custom canine-array 

We designed a custom array containing an average 
probe spacing of 3.3 Kb across the whole canine genome 
for CNV screening. Before analyzing potential CNVs 
among the seven cloned dogs, we first validated 
whether this custom array platform could reliably de- 
tect CNVs between different canine breeds. For this, we 
screened CNVs between the nuclear donor, a Labrador 
retriever genome (male), and a Boxer genome as refer- 
ence (female) using array-CGH, and verify the CNVs by 
genomic qPCR. Under our CNV detection criteria, we 
identified 72 CNVs between the Labrador retriever and 
Boxer across diverse chromosomes spanning 0.35% 
(8.6 Mb) of the canine genome. The average CNV size 
was 119 Kb (ranging from 11.1 to 450.6 Kb). Figure 1A 
illustrates the whole-genome plot of copy number gain- 
and loss-CNVs (CNV-G and CNV-L, respectively) identi- 
fied with our array. Details of the 72 CNVs are available in 
Additional file 1: Table SI. When we compared the 72 
CNVs with previously identified dog CNV regions [14-17], 
90.3% (65 of 72 CNVs) identified in this study overlapped 
with the previously reported canine CNVs, while seven 
(9.7%) were novel. Statistical significance of enriched CNV 
overlap with known dog CNV data was evaluated by per- 
mutation test as described previously [18]. The overlap 
with known dog CNV data was significantly not a random 
event (P < 0.001). To validate our array platform and our 
CNV defining criteria through genomic qPCR, we ran- 
domly selected nine CNV loci (2 CNV-Ls, and 7 CNV- 
Gs) out of the 72 CNVs and designed eleven primer sets 
targeting the nine loci. All the eleven qPCRs in the nine 
CNV loci showed consistent copy number changes that 
matched the array results (Figures IB, 1C). We also vali- 
dated the seven novel CNVs. The average size of the novel 



CNVs was 103.5 Kb (ranging from 26.8 to 235.7 Kb) and 
all of them were consistently validated by genomic qPCR 
(Additional file 1: Figure SI). Since one of the 7 novel 
CNVs (on chromosome 4) was included in the nine valid- 
ation targets, a total of 15 CNVs were qPCR validated in 
this study. All these data supported the reliability of our 
customized array platform and the subsequent analyses 
for identifying canine CNVs, therefore we applied this sys- 
tem to analyze cloned dog genomes. 

Association with segmental duplications and repetitive 
sequences 

To understand the mechanism of CNV formation, we 
assessed the association of CNVs with segmental dupli- 
cations (SDs) and repetitive sequences. When we com- 
pared the 72 CNVs with the CNV data from Nicholas et 
als report, which systematically analyzed SDs and associ- 
ated CNVs in canis familiaris [15], 48 (67%) overlapped 
with the SDs. Total length of the overlapped CNVs was 
about 6 Mb (-70% of the total length of the 72 CNVs). 
This result was coherent with Nicholas et als previous 
report demonstrating that 70% of CNVs were associated 
with SDs [16]. Considering that SDs are often substrates 
of CNV formation via non allelelic homologous recom- 
bination (NAHR) [3], ours and previous results suggest 
that SD-associated NAHR might be the major mechan- 
ism of canine CNV generation. When we observed the 
genomic fraction (fraction of retroelements per kb) of long 
interspersed elements (LINE), short interspersed elements 
(SINE) and long terminal repeats (LTR) within (intraCNV) 
or in the vicinity of the CNVs (up to 10 Mb), LINE ele- 
ments were relatively enriched the flanking regions (up to 
10 kb away from CNVs), while relatively depleted LINEs/ 
SINEs in intraCNV (Additional file 1: Figure S2). This re- 
sult is coherent with the previous reports suggesting that 
repeated sequences may play a role in the formation of 
structural variation and genomic diversity [19,20]. 

De novo CNVs identified in the cloned dogs 

To verify the genetic homogeneity between the donor 
dog and clones, we first performed the microsatellite - 
based evaluation of the individual identity. Being consist- 
ent with our previous report [12], all seven cloned dogs 
showed identical microsatellite patterns with the donor 
at all nine microsatellite markers examined (Additional 
file 1: Table S2). We calculated the probability that the 
clone should have the same genotype as the donor under 
two different assumptions: assuming Hardy- Weinberg 
equilibrium, and assuming that both individuals are sib- 
lings as described previously [21,22]. The probabilities of 
exact allele matching at the nine microsatellite markers 
in a population of Labrador retriever (Hardy- Weinberg 
P = 0.00048 and Full siblings P = 0.01763) suggest that 
allele matching at all nine microsatellite loci is not a 



Jung et al. BMC Genomics 2013, 14:863 
http://www.biomedcentral.com/1471 -21 64/1 4/863 



Page 3 of 1 1 



Mr 



**4 



I 5 I 6 I 7 I 



I 12 I 13 I 14 I 15 I 16 I 17 I 18 I 19 I 20 I 21 I 22 I 23 I 24 I 25 I 26 I 27 I 28 I 29 I 30 I 31 I 32 1 33 1 34 |35 1 36 1 37 |38 ' 



m 



0m 



B 



Chromosome 9 



Lh Pi L imifTi^iNji 



Chromosome 18 



-1 



-2t 



Chromosome 4 



18.1 19.7 
Chromosome 5 



22.9 32.1 32.8 
« Chromosome 15 



.1 



34.3 1.7 

2 
1 



Chromosome 19 

1 


: 




'23P25 


F 





83.2 32.2 



33.6 21.5 



Chromosome 23 



Chromosome 25 



Chromosome 32 



25.2 47.0 



49.1 40.2 



2.5- 
2.0- 
1.5 
1.0- 
0.5 
0.0- 



■ ■ 



I 



1 1 1 1 1 1 1 



P8 P9 P19 P6 P7 P14 P23 P25 P29 P30 P31 

Figure 1 Whole-genome plot of CNVs and qPCR validation. (A) Whole-genome plot of CNVs identified between Labrador retriever and Boxer 
genomes. Nine CNV regions (7 CNV-Gs, red arrow; 2 CNV-Ls, green arrow) were randomly selected for qPCR validation. The X-axis represents indi- 
vidual chromosomes and the Y-axis represents signal intensity ratios (Labrador retriever/Boxer) on a log2 scale. (B) Log2ratio plot around the nine 
selected CNV regions. Black bars represent primer position for qPCR validation. The X-axis represents genomic position (Mb) and the Y-axis repre- 
sents signal intensity ratios (Labrador retriever/Boxer) on log2 scale. (C) Three primer sets located in the CNV-L regions detected copy loss, and 
the other eight primer sets located in the CNV-G regions detected copy gain. The red and green line represent detection criterion for copy gain 
and copy loss, respectively. The Y-axis represents fold change based on the Boxer as the calibrator. 



random event. We next performed genome-wide array- 
CGH-based CNV screening to identify de novo CNVs in 
cloned dogs. For this, genomic DNA from the nuclear 
donor dog was used as the universal reference for each 
array- CGH. Therefore, all the CNVs identified can be 
interpreted as de novo ones. Under the same CNV 



detection criteria as in the validation experiment, we 
identified five de novo CNVs in two of the seven cloned 
dogs ranging from 36 Kb to 5.2 Mb (Table 1). Among 
the five de novo CNVs, a 36 Kb-sized CNV on chromo- 
some 2 (cloned dog-6) did not overlap with the previ- 
ously reported canine CNVs, while the other four 
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Table 1 De novo CNVs in cloned dogs 


Clone 


Chr 


Start (bp) 


End (bp) 


Length (bp) 


Event 


Genes* 


Validated by + 


Reported CNV § 


Clone-3 


15 


32,983,895 


33,026,044 


42,150 


CN Gain 




Q,S 


Yes 


Clone-3 


19 


23,000,913 


23,338,557 


337,644 


CN Gain 




Q, S, B 


Yes 


Clone-3 


32 


41,654,990 


41,731,424 


76,435 


CN Gain 




Q,S 


Yes 


Clone-3 


X 


1,401,988 


6,598,355 


5,196,367 


CN Gain 


ARSE, ARSH, OBP, GYG2, ARSD, ARSF, MXRA5, 
PRKX, PRKX NLGN4T, NLGN4X, HDHD 1, STS, 
PNPLA4, AKL1, TAF9, APOO, SHROOM2 


Q, S, B 


Yes 


Clone-6 


2 


71,652,481 


71,688,546 


36,066 


CN Loss 


ZBTB8A, PCNAP1 


Q 


No 



*Bold text indicates reference genes in canine and regular text indicates reference genes in human. 
tQ, genomic qPCR; S, dye-swap array-CGH analysis; B, B-allele profile analysis. 
§Overlapped with the previously reported CNVs [14-17], 



(cloned dog-3) overlapped with those reported (Table 1). 
The overlap with known dog CNV data was significantly 
not a random event (P < 0.001). Figure 2 A is the genome- 
wide signal intensity plot (log2 scale) in cloned dog-3 
(clone-3) as an example. To validate the de novo CNVs, 
we first carried out dye-swap array-CGH analysis for the 
clone-3 genome, but not for the clone-6 genome due to 
the lack of genomic DNA. Dye swap is a repeating 
hybridization on two-dye microarrays with the same 



samples but with swapped fluorescent labels. If there is 
true copy number change, the test/ reference intensity ra- 
tio value must be inverted by the switching of the dyes. 
The concordant copy number changes detected in dye- 
swap analysis were considered to be true changes. As ex- 
pected, all four CNVs identified in clone-3 consistently 
showed flipped signal intensity plots in the dye-swap 
hybridization (Figure 2B), supporting the accuracy of the 
CNVs identified by array-CGH. 



***** 
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Figure 2 Dye-swap analysis for the clone-3 genome. (A) Genome-wide signal intensity plot from clone-3. Four CNV regions were identified in 
the clone-3 genome (red arrow). The X-axis represents individual chromosomes and the Y-axis represents signal intensity ratio (clone-3/donor) on a 
log2 scale. (B) Log2ratio plots around the four CNV regions in clone-3. All CNV-Gs identified in clone-3 (red line) show flipped signal intensity plots 
(green line) in the dye-swap hybridization. The X-axis represents genomic position (Mb), and the Y-axis represents signal intensity ratio on a log2 scale. 
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B-allele profile of the de novo CNVs 

In addition to the dye-swap analysis, we also performed 
whole-genome SNP genotyping for the genomes of the 
donor and cloned dogs using an Illumina CanineHD 
170 K SNP array to obtain a more reliable interpretation 
of the de novo CNVs. A 5.2 Mb-sized region of copy 
number gain on chromosome X identified by array-CGH 
in clone-3 clearly showed two heterozygous clusters of 
SNPs in the B-allele plot, while no copy number gains or 
heterozygous clusters of SNPs were detected in the 
donor genome (Figure 3A). The 338 Kb-sized copy 
number gain on chromosome 19 of clone-3 also 
showed two heterozygous clusters of SNPs in the B- 
allele plot that were not detected in the donor genome 
(Figure 3B). The other three CNVs on chromosome 2, 
15 and 32 were too small to be interpreted by B-allele 
pattern (Additional file 1: Figure S3). 

In addition to the de novo CNVs, we identified a copy 
neutral de novo loss of heterozygosity (LOH) event on 
chromosome 2 by the B-allele pattern analysis. This 



37 Mb-sized LOH was identified in the distal end of 
chromosome 2 in clone-6, but not in the nuclear donor 
nor in the other cloned dogs (Figure 4). 

Genomic qPCR validation of the de novo CNVs 

To validate the existence and boundaries of the CNVs, we 
performed genomic qPCR for multiple points within and 
outside the expected CNV breakpoints. For this, we de- 
signed two to seven target specific primer pairs for each 
CNV (Additional file 1: Table S3). Through multiple-point 
genomic qPCR, four of the five CNVs were consistently 
validated, while the other one relatively smaller sized 
CNV was not validated (Figure 5). In the 5.2 Mb-sized 
region of copy number gain on chromosome X and the 
338 Kb-sized gain on chromosome 19 in clone-3, all 
the qPCRs within the CNV regions consistently showed 
copy number gains and the qPCRs outside CNV region 
showed diploid copy numbers, which confirmed the ac- 
curacy of our array-CGH analysis in defining boundar- 
ies (Figures 5 A, 5B). Likewise, all the other cloned dogs 



o 


1 


2 




CD 


0 








-1 




-2 




1.0 






c 

CD 
=5 


0.5 


CT 




2 

LL 


0.0 


CD 


1.0 


CD 


< 
CD 


0.5 




00 



Chromosome X 



126 (Mb) 




B 



0 


6 


19 (Mb) 




Chromosome 19 






— 

















18 



33 



56 (Mb) 




33 (Mb) 

Figure 3 Examples of allelic imbalance in clone-3. (A) Upper panel, Log2ratio plot on chromosome X in clone 3. A 5.2-Mb copy number gain 
region on chromosome X was identified in clone-3 (black arrow). The X-axis represents genomic position (Mb) and the Y-axis represents signal in- 
tensity ratio on a log2 scale. Lower panel, B-allele plot of the proximal part of chromosome X (dashed box of upper panel). A heterozygous cluster 
was detected in chromosome X (black arrow), and this allelic imbalance has the same genomic position as the copy gain region. (B) Upper panel, 
Log2ratio plot on chromosome 19 in the clone-3. A 338-Kb copy number gain region on chromosome 19 was identified in clone-3 (black arrow). 
Lower panel, B-allele plot of the dashed box on chromosome 19. A heterozygous cluster was detected on chromosome 19 (black arrow), and this 
allelic imbalance has same position as the copy gain region. 
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Figure 4 Example of the LOH event in clone-6. A 37-Mb de novo LOH event was identified on the distal end of chromosome 2 in clone-6 that 
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containing no CNVs were diploid at all qPCR points. 
Of the relatively smaller-sized CNVs, a 36 Kb copy 
number loss on chromosome 2 in clone-6 and 42 Kb 
gain on chromosome 15 in clone-3 were successfully 
validated and their boundaries were estimated by qPCR 
(Figure 5C, 5D). However, the 76 Kb-sized CNV-G on 
chromosome 32 in clone-3 was not validated by gen- 
omic qPCR (Figure 5E). 

Discussion 

In this study, we aimed to identify de novo post-cloning 
CNV events which may exist in cloned dogs like in hu- 
man MZ twins. For this purpose, we designed a 720 K 
custom array by merging the NimbleGen 385 K dog 
array as a backbone and 335 K probes for known canine 
CNV regions [14]. Using this platform, we identified five 
de novo CNVs in two of the seven cloned dogs from the 
same nuclear donor. Our findings support the occur- 
rence of de novo post-cloning CNV events in mammals 



including dogs. Although de novo CNV events in 
phenotypically concordant MZ twins have been reported 
in humans [8], to the best of our knowledge, this is the 
first report of de novo CNVs in cloned animals including 
dogs generated by SCNT technology. 

As mentioned above, the de novo CNVs were unevenly 
identified across the dog clones (2/7 clones); four CNVs 
were identified in one clone (clone-3) and one was in 
another clone (clone-6). Our result is similar with the 
previous reports observing the CNVs in human MZ 
twins. Breckpot et al. identified CNVs in 1 out of 6 
phenotypically discordant MZ twins [9]. In Ehli et al's 
report with 50 concordant and discordant MZ twins, 18 
de novo post-twining CNVs were detected in 14 out of 
45 twin pairs [23]. However, we cannot rule out the pos- 
sibility of existence of more de novo CNVs because the 
array platform and defining algorithm we used in this 
study cannot be perfect. Also, it is possible that some 
genetic events may have happened in the clone-3 
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harboring 4 CNVs during its development or SCNT 
procedure. Further CNV studies with a larger number 
of cloned dogs will be required to understand this 
phenomenon. 

When we screened CNVs between the donor Labrador 
retriever and a Boxer to validate the performance of our 



custom array, the number of CNVs was largely similar 
to the report by Nicholas et al. that compared CNVs 
using a 2.1 M array between nine breeds of modern do- 
mestic dogs including Labrador retriever [16]. Most of 
the CNVs (90.3%) identified in our validation experi- 
ment overlapped with the previously reported CNVs 
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[14-17], and the overlap was significantly not random 
event (P < 0.001). All the nine randomly selected CNV 
loci and seven novel CNVs identified by array-CGH 
showed consistent copy number changes that were vali- 
dated by qPCR. These results suggest the reliability of our 
custom array and CNV defining procedure. Although ex- 
ploring the novel canine CNVs is not our aim in this 
study, this result suggests that a substantial number of ca- 
nine CNVs that have yet to be identified. 

To validate the five de novo CNVs in the cloned dogs, 
we performed three different verification experiments; 
genomic qPCR, dye-swap array-CGH analysis and B- 
allele profile analysis. Two larger de novo CNVs (5.2 Mb 
and 338 Kb) on chromosomes X and 19 in clone-3 were 
consistently validated by all three experiments. The 
other three CNVs on chromosomes 2, 15, and 32 (36.1 - 
76.4 Kb) seem to be too small to be interpreted by B-allele 
pattern (Additional file 1: Figure S3). According to the 
technical note about CNV analysis algorithm for Illumina 
HumanHap240S array, which has a higher SNP density 
than the array used in this study (170 K), the minimum 
recommended window size for LOH score algorithm is 
360 Kb [24]. Of the three smaller-sized CNVs, two (a CNV 
on chromosome 2 in clone-6 and another one on chromo- 
some 15 in clone-3) were verified by qPCR, but the CNV 
on the distal end of chromosome 32 was not verified. 
However, the CNV region on chromosome 32 was vali- 
dated in dye-swap analysis, which suggests that the 76 Kb- 
sized CNV region on chromosome 32 might be real. 

There are several genes in the 5.2 Mb-sized copy 
number gain CNV locus including the neuroligin 4 X- 
linked (NLGN4X), steroid sulfatase isozyme S (STS), 
and odorant-binding protein (OBP) genes. In human, 
genetic variants of the NLGN4X gene have been re- 
ported to be associated with autism spectrum disorder 
in Japanese and Han Chinese populations [25-27]. A 
CNV of the STS gene (duplication) has also been re- 
ported to be associated with intellectual disability in 
human [28]. OBP is small soluble carrier protein that 
belongs to the family of lipocalins and that may be in- 
volved in the perception of pheromones [29]. However, 
we did not observe any differences regarding the phe- 
notypes related to these genes in the cloned dogs har- 
boring the de novo CNVs in this study. 

In addition to the de novo CNVs, we also detected a 
copy neutral de novo LOH event on chromosome 2. 
This 37 Mb-sized LOH event occurred in a cloned dog 
but was not present in the nuclear donor dog, indicating 
that the uniparental disomy (UPD) arose from a mitotic 
event in the cloned dog. It is well known that UPD is a 
genetic event associated with human diseases and mosai- 
cism of aneuploidy can be associated with UPD [30]. 
Our results suggest that SCNT could result in UPD- 
related phenotypic differences among clones. 



There are several limitations in this study. First, since 
the sample size of this study is very small, the scope of 
our analysis could not go beyond describing de novo 
CNVs. Second, even though we identified de novo CNVs 
in cloned dogs, we could not analyze any association be- 
tween the CNVs and phenotypic characteristics, because 
there have been no phenotypic differences among them 
so far. To overcome both limitations and draw biologic- 
ally meaningful conclusions, we need to clone more dogs 
and to carefully observe the phenotypes of them, and 
apply more informative technologies such as whole gen- 
ome sequencing. 

Conclusions 

In conclusion, we identified de novo CNV and LOH 
events in cloned dogs generated from one nuclear 
donor. To the best of our knowledge, this is the first re- 
port of de novo CNVs in cloned animals including dogs 
generated by SCNT technology. Since cloned animals 
share almost identical genetic background like MZ 
twins, to study their de novo genetic events can help 
understand formation mechanisms of genetic variants 
and their biological implications. Also, our results can 
be a useful tool to study the genetic basis of diverse dis- 
ease/phenotype traits in dogs. 

Methods 

Source of dog DNA samples 

Seven cloned Labrador retriever dogs named clone- 1 
to -7 were generated by somatic cell nuclear transfer 
(SCNT) using fibroblasts derived from a seven-year- 
old Labrador retriever, and the genetic homogeneity was 
validated by microsatellite analysis using nine markers 
(PEZ01, PEZ02, PEZ03, PEZ06, PEZ13, PEZ17, FH2079, 
FH2054, FH2010) in the previous study [12]. In order to 
extract genomic DNA samples, blood was collected from 
the seven cloned dogs at three years of age and from the 
donor dog at ten years of age. Approximately 5 ml of 
blood was collected into anticoagulant tubes containing 
EDTA. Genomic DNA was extracted from whole-blood 
using a DNeasy Blood & Tissue Kit (Qiagen, Hilden, 
Germany) and quantified by NanoDrop spectrophotom- 
eter. This study was performed with the approval of ani- 
mal ethics committee of Seoul National University. 

Whole-genome array-CGH experiment 

We used a NimbleGen custom 720K canine whole- 
genome array platform (Roche NimbleGen, Penzberg, 
Germany), designed from the UCSC dog genome build 
of May 2005 (Broad/canFam2). This array is comprised 
of 385,000 probes evenly distributed throughout the 
unique sequence of the genome as a backbone and an 
additional -330,000 probes which were designed based 
on previously reported dog CNV regions [14]. The 
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average probe spacing for our custom array platform 
was -3.3 Kb. The array-CGH experiments were per- 
formed according to the manufacturer s instructions. In 
brief, 1 \ig genomic DNA from each cloned dog was la- 
beled with Cy3-dCTP. The reference DNA was labeled 
with Cy5-dCTP. Each dye-labeled DNA was purified and 
mixed with hybridization reagents, applied on the array 
and incubated for 48 hours at 42°C in a MAUI 
hybridization machine (BioMicro, Salt Lake City, UT). 
After washing the slides, arrays were scanned using a 
GenePix 4000B scanner (Axon Instruments, Sunnyvale, 
CA). Image separation, signal intensity extraction, and 
normalization were performed using the NimbleScan 
software version 2.5 (Roche NimbleGen). 

Detection of CNVs 

The Rank-Segmentation CNV calling algorithm in 
NEXUS software v3.1 (BioDiscovery Inc., El Segundo, 
CA) was used to define CNVs of each sample. The pa- 
rameters for defining CNVs were as follows: significance 
threshold = 1.0E-6; maximum contiguous probe spacing 
(Kbp) = 1000 Kb; minimum number of contiguous 
probes per CNV segment = 5; threshold of signal inten- 
sity ratio > 0.2 on the log2 scale for gains and < -0.2 on 
the log2 scale for losses. 

Whole-genome SNP genotyping 

Genome-wide SNP genotyping was conducted using the 
CanineHD SNP array, which contains 173,662 SNP 
markers (Illumina, Inc., San Diego, CA) according to 
manufacturers instructions. The average probe spacing 
is 13 Kb. The raw data consisted of the signal intensity 
(log R ratio: LRR) and allelic intensity (B allele fre- 
quency) that were obtained using GenomeStudio soft- 
ware (Illumina, Inc.). Based on the B allele frequency 
profile, we defined the CNVs according to the concept of 
integrating SNP and CNV genotyping [3]. Samples that 
did not have a standard deviation (SD) of LRR < 0.24 were 
excluded for quality control. The SNPRank Segmentation 
algorithm in NEXUS software was used to define allele 
imbalance of each sample. 

Genomic quantitative PCR analysis 

Array-CGH results were experimentally validated by 
genomic qPCR in randomly selected CNV loci for boxer 
and all CNV loci for cloned dogs. As a diploid internal 
control, a genomic region on chr7: 18946362 that showed 
no genomic alteration in the array-CGH data was used. 
Details including primer information for targets and the 
diploid control locus are available in Additional file 1: Table 
S3. Genomic qPCR was performed using the ViiA7™ real- 
time PCR system (Life Technologies, Carlsbad, CA) as de- 
scribe elsewhere [31]. In brief, a 10 \i\ of reaction mixture 
contained 10 ng of genomic DNA, THUNDERBIRD™ 



SYBR qPCR Mix (TOYOBO, Osaka, Japan), lx ROX, and 
6 pmole of each primer. Thermal cycling conditions con- 
sisted of one cycle of 1 min at 95°C followed by 40 cycles 
of 5 sec at 95°C, 10 sec at 61°C, and 20 sec at 72°C. All 
qPCR experiments were triplicate and the copy number of 
each target was defined as 2" AACT , where AC t is the differ- 
ence in threshold cycles of the sample in question normal- 
ized against the reference region and expressed relative to 
the value obtained by calibrator DNA (individual/calibrator) 
as describe elsewhere [32]. DNA from the nuclear donor 
dog was used as calibrator for all qPCR experiment. 

Probability of identity 

Probability of identity (PI) for multilocus microsatellite 
data was computed using GENECAP program [33] 
under the two different assumptions of Hardy- Weinberg 
equilibrium and the sibling PI as described previously 
[21,22]. Where, PI indicates the probability of two indi- 
viduals within the population sharing the same genotype. 
In GENECAP, multilocus PI values were computed by 
mutiplying all locus-specific Pis under the assumption 
that each locus is independent. 

Significance of enriched CNV overlap with known dog 
CNV data 

Statistical significance of enriched CNV overlap with 
known dog CNV data was evaluated by permutation 
test as described previously [18]. In this test, original 
CNV overlap with database was calculated based on >1 bp 
overlap measure. In each iteration, CNVs were ran- 
domly reshuffled while maintaining their sizes and 
numbers within each chromosome. Based on these per- 
muted CNVs, we calculated the overlap with the same 
database. After repeating this process 1000 times, we 
calculated P-value as the fraction of tests where the 
number of overlaps was greater than the original unper- 
muted CNV overlap. 

Statistical analyses 

Statistical analyses were performed using SPSS version 
18.0 (SPSS Inc., Chicago, IL). P values < 0.05 were con- 
sidered significant. 

Availability of supporting data 

The array-CGH data (NimbleGen custom 720 K canine 
whole-genome array) and the SNP genotyping data 
(Illumina CanineHD SNP array) from this study are 
submitted to the NCBI Gene Expression omnibus 
(http://www.ncbi.nlm.nih.gov/geo) under accession no 
GSE49092 and GSE49123, respectively. The following 
supporting data are available with the online version of 
this paper. 
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Additional file 



Additional file 1: Table SI. CNVs between Labrador retriever and Boxer 
genome. Table S2. Microsatellite analysis of the donor cell and cloned dogs 
using nine canine specific microsatellite markers. Table S3. Sequence 
information of the primers for genomic qPCR validation. Figure SI. Novel 
CNVs identified between a Labrador retriever and a Boxer, and qPCR 
validation. Figure S2. Genomic fraction of the repetitive sequence elements. 
Figure S3. B allele profiles around the three small-sized de novo CNVs. 
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